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Abstract 

The reductive perturbation method has been employed to derive the 
Korteweg-de Vries (KdV) equation for small but finite amplitude elec- 
trostatic waves. The Lagrangian of the time fractional KdV equation is 
used in similar form to the Lagrangian of the regular KdV equation. The 
variation of the functional of this Lagrangian leads to the Euler-Lagrange 
equation that leads to the time fractional KdV equation. The Riemann- 
Liouvulle definition of the fractional derivative is used to describe the 
time fractional operator in the fractional KdV equation. The variational- 
iteration method given by He is used to solve the derived time fractional 
KdV equation. The calculations of the solution with initial condition 
Ao sec h(cx) 2 are carried out. Numerical studies have been made using 
plasma parameters close to those values corresponding to the dayside au- 
roral zone. The effects of the time fractional parameter on the electrostatic 
solitary structures are presented. 

Keywords: Euler-Lagrange equation, Riemann-Liouvulle fractional deriva- 
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1 Introduction 

Because most classical processes observed in the physical world are nonconser- 
vative, it is important to be able to apply the power of variational methods 
to such cases. A method used a Lagrangian that leads to an Euler-Lagrange 
equation that is, in some sense, equivalent to the desired equation of motion. 
Hamilton's equations are derived from the Lagrangian and are equivalent to the 
Euler-Lagrange equation. If a Lagrangian is constructed using noninteger-order 
derivatives, then the resulting equation of motion can be nonconservative. It 
was shown that such fractional derivatives in the Lagrangian describe noncon- 
servative forces [1, 2]. Further study of the fractional Euler-Lagrange can be 
found in the work of Agrawal [3, 4], Baleanu and coworkers [5, 6] and Tarasov 
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and Zaslavsky [7, 8]. During the last decades, Fractional Calculus has been 
applied to almost every field of science, engineering and mathematics. Some 
of the areas where Fractional Calculus has been applied include viscoelasticity 
and rheology, electrical engineering, electrochemistry, biology, biophysics and 
bioengineering, signal and image processing, mechanics, mechatronics, physics, 
and control theory [9] . 

On the other hand, electron acoustic waves (EAWs) propagation in plasmas 
has received a great deal attention because of its vital role in understanding 
different types of collective processes in laboratory devices [10, 11] as well as 
in space environments [12, 13]. It has been argued that when the hot to cold 
electron temperature ratio is greater than 10, the electron-acoustic mode may 
be the principal mode of the plasma in which the restoring force comes from the 
pressure of the hot electrons, while the inertia comes from the mass of the cold 
electron component [14]. The ions play the role of a neutralizing background, 
i.e., the ion dynamics does not influence the EAWs because its frequency is 
much larger than the ion plasma frequency. Several theoretical attempts have 
been made to explain nonlinear EAWs in plasma systems [15-17]. 

To the author's knowledge, the problem of time fractional KdV equation in 
collisionless plasma has not been addressed in the literature before. So, our mo- 
tive here is to study the effects of time fractional parameter on the electrostatic 
structures for a system of collisionless plasma consisting of a cold electron fluid 
and non-thermal hot electrons obeying a non-thermal distribution and station- 
ary ions. Our choice of non-thermal distribution of electrons is prompted by 
its convenience rather than as precise fitting of the observations. We expect 
that the inclusion of the non-thermal electrons, time fractional parameter will 
change the properties as well as the regime of existence of solitons. 

Several methods have been used to solve fractional differential equations 
such as: the Laplace transform method, the Fourier transform method, the 
iteration method and the operational method [18, 19]. Recently, there are some 
papers deal with the existence and multiplicity of solution of nonlinear fractional 
differential equation by the use of techniques of nonlinear analysis [20-23]. In this 
paper, the resultant fractional KdV equation will be solved using a variational- 
iteration method (VIM) firstly used by He [24, 25]. 

This paper is organized as follows: Section 2 is devoted to describe the 
formulation of the time-fractional KdV (FKdV) equation using the variational 
Eulcr-Lagrange method. In section 3, the resultant timc-FKdV equation is 
solved approximately using VIM. Section 4 contains the results of calculations 
and discussion of these results. 

2 Basic Equations and Time-fractional KdV equa- 
tion 

We consider a homogeneous system of unmagnctized collisionless plasma con- 
sisting of cold electrons fluid, non-thermal hot electrons obeying a non-thermal 
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distribution and stationary ions. Such system is governed by the following nor- 
malized equations in one dimension [26]: 

d d 

— u c (x, t) + —[u c (x, t) n c (x, t)} = 0, (la) 

Q Q 

— u c (x, t)+u c (x, t)—u c (x, t)-jd(x, i) = 0, (lb) 
with Poisson's equation 

d 2 1 1 

— (j>(x, t) - -n c (x, t) - n h (x, t) + (1 + -) = 0. (lc) 
ox* 7 7 

The non-thermal hot electrons density n^(x, t) is given by: 



n h (x, t) = [l-/3 <f>(x, t)+0 6{x, t) 2 ] ex P [0(x, t)\, = 45/(1 + 3J). (2) 

In these equations, n c (x, t)[nh(x, t)] is the cold (non-thermal hot) electrons 
density normalized by equilibrium value n c o[nho], u c (x, t) is the cold electrons 

fluid velocity normalized by hot electron acoustic speed C e = , 4>{x, t) 

is the electric potential normalized by ka J H , 7 = is the hot to cold electron 
equilibrium densities ratio, m e is the electron mass, <5 is a parameter which de- 
termines the population of energetic non-thermal hot electrons, e is the electron 
charge, x is the space co-ordinate and t is the time variable. The distance is 
normalized to the hot electron Debye length Xoh, the time is normalized by the 
inverse of the cold electron plasma frequency cj^ 1 and ks is the Boltzmann's 
constant. Equations (la) and (lb) represent the inertia of cold electron and (lc) 
is the Poisson's equation needed to make the self consistent. The hot electrons 
are described by non-thermal distribution given by (2). 

According to the general method of reductive perturbation theory, we intro- 
duce the slow stretched co-ordinates [27]: 

T = e 2 / 3 t and£ = e 1/2 (a;-At), (3) 

where e is a small dimensionless expansion parameter and A is the wave speed 
normalized by C e . All physical quantities appearing in (1) are expanded as 
power series in e about their equilibrium values as: 

n c (f , t) = 1 + eni(f , r) + e 2 n 2 (£, r) + e 3 n 3 (£, r) + (4a) 
u c (£, r) = eui(£, r) + e 2 u 2 (^ r) + e 3 u 3 (£, t) + (4b) 

</>(£, r) = efafo t) + e 2 2 (£, r) + e 3 ^, r) + .... (4c) 

We impose the boundary conditions as £ — > 00, n c = rih = 1, u c = and 
6 = 0. 
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Substituting (3) and (4) into (1) and equating coefficients of like powers of 
e, the lowest-order equations in e lead to the following results: 



m(£,T) = ±MZ,t) and Ui (£,t) = {MS,t). (5) 



-^M^t) and u i(£, T ) = J 
Poisson's equation gives the linear dispersion relation 



Considering the coefficients of 0(e 2 ) and eliminating the second order-perturbed 
quantities n 2 , u 2 and </> 2 lead to the following KdV equation for the first-order 
perturbed potential: 



^MZ, r) + AM^ r)-^M^ r) + B |^<M(. r) = 0. (7a) 

where 



<fii(£, t) is a field variable, £ is a space coordinate in the propagation direction 
of the field and r G T(= [0, To]) is the time. The resultant KdV equation (7a) 
can be converted into time-fractional KdV equation as follows: 

Using a potential function v(£, r) where <j>i(£, r) = v^(^, t) gives the 
potential equation of the regular KdV equation (1) in the form 

vtrfo t) + A 7>«(£, r) + B vkk(S, r) - 0, (8) 

where the subscripts denote the partial differentiation of the function with re- 
spect to the parameter. The Lagrangian of this regular KdV equation (7) can 
be defined using the semi-inverse method [28, 29] as follows. 

The functional of the potential equation (8) can be represented by 



J(v) = JdZj d T {v(£, T )[civt T (£,T) + c 2 Avs(Z,t)vk(£,t) + c 3 5^(£, r)]}, 

R T 

(9) 

where ci, c 2 and C3 are constants to be determined. Integrating by parts and 
taking v t \r = v^\r — v^\t — lead to 



R T 



(10) 

The unknown constants c, (i — 1, 2, 3) can be determined by taking the 
variation of the functional (10) to make it optimal. Taking the variation of this 
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functional, integrating each term by parts and make the variation optimum give 
the following relation 

2ci« er (£,T) +3c 2 Av^,t)v^,t) +2c 3 Bv m ^,r) = 0. (11) 

As this equation must be equal to equation (8), the unknown constants are 
given as 

ci = 1/2, c 2 = 1/3 and c 3 = 1/2. (12) 

Therefore, the functional given by (10) gives the Lagrangian of the regular 
KdV equation as 

L(v T , vt,vx) = ~v^,t)v t ^,t) - IavK^t) + ^Bv 2 ^,t). (13) 

Similar to this form, the Lagrangian of the time-fractional version of the 
KdV equation can be written in the form 

F( D°v, v vx) = -l[oDy^,T)}v^,T)-^Avl^,r) + ^Bvl^,r), 

< a<\, (14) 

where the fractional derivative is represented, using the left Riemann-Liouville 
fractional derivative definition as [18, 19] 



a B?f(t) = T ^_ a J &\f dTit-r)"-"- 1 /^)], k-l < a < 1, t e [a, b}. (15) 
The functional of the time-FKdV equation can be represented in the form 
J(v) = fdtij dTF( D«v, w 6 w«), (16) 

R T 

where the time-fractional Lagrangian F(oD^v, v^,v^) is defined by (14). 

Following Agrawal's method [3, 4], the variation of functional (16) with 
respect to r) leads to 

SJ(v) = / dr{^- v 5 oD> + + °Lsv«}. (17) 

R T 

The formula for fractional integration by parts reads [3, 18, 19] 

f dtf(t) a D?g(t)= f dtg(t) t DZf(t), f(t),g(t)e[a,b}. (18) 

J a J a 

where t D^ , the right Riemann-Liouville fractional derivative, is defined by [18, 
19] 
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tD%f{t) = AJ¥- *.[J dr(T-t) k - a - 1 f{r)], fc-l<«<M£ [a,b]. (19) 
Integrating the right-hand side of (17) by parts using formula (18) leads to 

^y^Wse^^spW 14 "' <20) 

R T 

where it is assumed that 5v\t — Sv\r = 5v^\r = 0. 

Optimizing this variation of the functional J(v), i. e; SJ(v) — 0, gives the 
Euler-Lagrange equation for the time-FKdV equation in the form 

^<5g;> < 21 » 

Substituting the Lagrangian of the time-FKdV equation (14) into this Euler- 
Lagrange formula (21) gives 



- 2 TDZ vt& r) + - D?v&, r) + Av t {£, r)v t ^, r) + Bv mi {^ r) = 0. (22) 

Substituting for the potential function, v^(^,t) — 4>i(£,t) — $(£,r), gives 
the time-FKdV equation for the state function $(£, r) in the form 



-[ D?$(Z,t) - T r)] + A<D(e,r)$ e (e,r) + 5$ ecc (£,r) = 0, (23) 

where the fractional derivatives n D" and T D^ g are, respectively the left and 
right Riemann-Liouville fractional derivatives and are defined by (15) and (19). 
The time-FKdV equation represented in (14) can be rewritten by the formula 

\ R D?$(£, r)+A $(£, r)$ e (£, r) + B $ eec (£, r) = 0, (24) 

where the fractional operator qD" is called Riesz fractional derivative and can 
be represented by [4, 18, 19] 



oA°7(*) = \[oD?f{t)+ (-l) k t DZJ(t)} 

fc-l < a<l,t€[a,b]. (25) 

The nonlinear fractional differential equations have been solved using dif- 
ferent techniques [18-23]. In this paper, a variational-iteration method (VIM) 
[24, 25] has been used to solve the time-FKdV equation that formulated using 
Euler-Lagrange variational technique. 



6 



3 Variational-Iteration Method 



Variational-iteration method (VIM) [24, 25] has been used successfully to solve 
different types of integer nonlinear differential equations [30, 31]. Also, VIM is 
used to solve linear and nonlinear fractional differential equations [32, 33]. This 
VIM has been used in this paper to solve the formulated time-FKdV equation. 

A general Lagrange multiplier method is constructed to solve non-linear 
problems, which was first proposed to solve problems in quantum mechanics 
[24]. The VIM is a modification of this Lagrange multiplier method [25]. The 
basic features of the VIM are as follows. The solution of a linear mathematical 
problem or the initial (boundary) condition of the nonlinear problem is used as 
initial approximation or trail function. A more highly precise approximation 
can be obtained using iteration correction functional. Considering a nonlinear 

partial differential equation consists of a linear part LU(x,t), nonlinear part 

NU(x,t) and a free term f(x,t) represented as 

LU(x,t) + NU(x,t) = f(x,t), (26) 

where L is the linear operator and N is the nonlinear operator. According to the 
VIM, the (n + l)th approximation solution of (26) can be given by the iteration 
correction functional as [24, 25] 

U n+1 (x,t) = U n (x,t)+ [ dT\{T)[LU n (x,T)+NU n (x,T)-f(x,T)},n>0, (27) 

where A(r) is a Lagrangian multiplier and U n (x,T) is considered as a restricted 
variation function, i. e; SU n (x, r) = 0. Extreme the variation of the correction 
functional (27) leads to the Lagrangian multiplier A(r). The initial iteration 
can be used as the solution of the linear part of (26) or the initial value U(x, 0). 
As n tends to infinity, the iteration leads to the exact solution of (26), i. e; 

U(x,t) = lim U n (x,t). (28) 

n— >oo 

For linear problems, the exact solution can be given using this method in 
only one step where its Lagrangian multiplier can be exactly identified. 

4 Time- fractional KdV equation Solution 

The time-FKdV equation represented by (24) can be solved using the VIM by 
the iteration correction functional (27) as follows: 

Affecting from left by the fractional operator on (24) leads to 
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-*D t 1 -°[A*&t)^*&t)+B^$(Z,t)], 
< a < 1, r G [0,T ], (29) 
where the following fractional derivative property is used [18, 19] 



fD 6 -[*D//(t)] = fDr"/(t)-x;?D/-V(t)u|rJ" J 

fc — 1 < j3 <k. (30) 



As a < 1, the Riesz fractional derivative ^D T a 1 is considered as Riesz 
fractional integral qI^" that is defined by [18, 19] 

*J T "/(t) = i[o/ T "fW + T I 6 °7(t)] = J" dr\t-rr\f(r), a>0. (31) 

where oI T a f(t) and T I b a f(t) are the left and right Riemann-Liouvulle fractional 
integrals, respectively [18, 19]. 

The iterative correction functional of equation (29) is given as 

*n+l(£,T) = $ n (£,T)+| T dT'A(T'){^$ n (£,T') 

r) ~ f) 3 ~ 

+ *D T )~ a [A $„(£,t') J> n (Z, T ') + B -3$nK,r')]},(32) 



where n > and the function t) is considered as a restricted variation 

function, i. e; <5<f> n (£,r) = 0. The extreme of the variation of (32) using the 
restricted variation function leads to 

<J$„+i(£,t) = 6$ n (Z,T)+ £ dr'\(T') ^$„K,t') 

= <5$„(£, r) + A(r) r) - £ dr'-^X(r') 5$ n (£, r') = 0. 

This relation leads to the stationary conditions 1+A(r) = and ^p- A(r') = 0, 
which leads to the Lagrangian multiplier as A(r') = —1. 
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Therefore, the correction functional (32) is given by the form 



*n+i(£,T) = $„(£,t)- f dr'{^7*nK,r') 

+ ^-"L4 * n (e,r')^*™(^r') + B^*nK,r')]},(33) 

where n > 0. 

In Physics, if r denotes the time-variable, the right Riemann-Liouville frac- 
tional derivative is interpreted as a future state of the process. For this reason, 
the right-derivative is usually neglected in applications, when the present state 
of the process does not depend on the results of the future development [3]. 
Therefore, the right-derivative is used equal to zero in the following calcula- 
tions. 

The zero order correction of the solution can be taken as the initial value of 
the state variable, which is taken in this case as 

< fo(C,T) = $(C,0) = ,4 sech 2 ( c e). (34) 
where Aq and c are constants. 

Substituting this zero order approximation into (33) and using the definition 
of the fractional derivative (25) lead to the first order approximation as 

$i(£,t) = AosechV) 

+2A c sinh(c£) sech 3 (c£)[4c 2 B 

+(A A- 12c 2 B) sec h 2 (c^)\ r /° ■ (35) 

1 (a + 1) 

Substituting this equation into (33), using the definition (25) and the Maple 
package lead to the second order approximation in the form 

* 2 (£,r) = A scch 2 (c£) 

+2A a c sinh(cC) sech 3 (c£)[4c 2 i? + (A A - 12c 2 B) sec /i 2 (c£)]— — ^ 

T(a + 1) 

(36) 

The higher order approximations can be calculated using the Maple or the 
Mathematica package to the appropriate order where the infinite approximation 
leads to the exact solution. 
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5 Results and calculations 



For small amplitude electron-acoustic waves, the time fractional Korteweg-de 
Vries equation has been derived. The Riemann-Liouvulle fractional derivative 
[18, 19] is used to describe the time fractional operator in the FKdV equation. 
He's variational-iteration method [24, 25] is used to solve the derived time-FKdV 
equation. 

To make our result physically relevant, numerical studies have been made 
using plasma parameters close to those values corresponding to the dayside 
auroral zone [15, 34]. 

However, since one of our motivations was to study effects of time frac- 
tional order (a), the energetic population parameter (S) and the hot to cold 
electron equilibrium densities ratio (7) on the existence of solitary waves. The 
electrostatic potential as a function of space variable £ and time variable r is 
represented in Figure (1) that has a single rarefactive soliton shape. In Figure 
(2), the relation between the hot to cold electron equilibrium densities ratio 
(7) and the amplitude of the electrostatic potential solitary wave $(0,r) is 
obtained at different values of the time variable r. It is seen that the soliton 
amplitude decreases with the increase of (7). In addition the soliton amplitude 
|$(0,t)| against the fractional order (a) is represented in Figure (3) at different 
time variable values. It is seen that |$(0,r)| increases with the increase of (a). 
Figure (4) shows the relation between |$(0,r)| and the energetic population 
parameter (S). The effect of the fractional order (a) on the electrostatic soliton 
amplitude |$(0,t)| for different values of the energetic population parameter 
(5) is given in Figure (5). To compare our result (the amplitude of the elec- 
trostatic potential |$(0,t)|) with that observed in the auroral zone, we choose 
a set of available parameters corresponding to the dayside auroral zone where 
an electric field amplitude E = 100 mV/m has been observed [15, 34] with 
T c = 5cV, Th — 250 eV, n c0 = 0.5 m -3 and riho = 2.5 m -3 . These parameters 
correspond to \oh ~ 7430 cm and the normalized electrostatic wave potential 
amplitude $0 = E k^r^ ~ 0.03 V, which is obtained for different values of (a) 
and (6) [cf. Figure (5)]. This value of the electrostatic potential amplitude is 
given at values of (a) and (5) as: (a — 0.61, 5 = 0.01), (a = 0.69, 5 = 0.05), 
(a = 0.78, S = 0.1) and (a = 0.9, S = 0.15) at velocity v = 0.04, 7 = 5 and time 
r = 10. 

In summery, it has been found that the amplitude of the electron acoustic 
solitary waves as well as the parametric regime where the solitons can exist are 
sensitive to the time fractional order. 

We have stressed out that it is necessary to include the space fractional 
parameter. This is beyond the scope of the present paper and it will be include 
in a further work in electron-acoustic solitary wave. The application of our 
model might be particularly interesting in the auroral region. 
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Figure Captions: 

Fig. (1): The electrostatic potential <&(£,t) against £ and r at 7 = 5, 
v = 0.04, a = 0.5 and 6 = 0.1. 

Fig. (2): The amplitude of the electrostatic potential |<I>(0,t)| against 7 at 
different values of time for v — 0.04, a = 0.5 and S = 0.1. 

Fig. (3): The amplitude of the electrostatic potential |$(0,r)| against a at 
different values of time for 7 = 5, v = 0.04 and 5 = 0.1. 

Fig. (4): The amplitude of the electrostatic potential |$(0,r)| against S at 
different values of time for 7 = 5, v = 0.04 and a — 0.5. 

Fig. (5): The amplitude of the electrostatic potential |$(0,r)| against a at 
different values of S for v = 0.04, r = 10 and 7 = 5. 
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